Nonequilibrium Steady State Driven by a Nonlinear Drift Force 
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We investigate the properties of the nonequilibrium steady state for the stochastic system driven 
by a nonlinear drift force and influenced by noises which are not identically and independently 
distributed. The nonequilibrium steady state (NESS) current results from a residual part of the 
drift force which is not cancelled by the diffusive action of noises. From our previous study for 
the linear drift force the NESS current was found to circulate on the equiprobability surface with 
the maximum at a stable fixed point of the drift force. For the nonlinear drift force, we use the 
perturbation theory with respect to the cubic and quartic coefficients of the drift force. We find an 
interesting potential landscape picture where the probability maximum shifts from the fixed point 
of the drift force and, furthermore, the NESS current has a nontrivial circulation which fiows ofl^ 
the equiprobability surface and has various centers not located at the probability maximum. The 
theoretical result is well confirmed by the computer simulation. 

PACS numbers: 87.23.Kg, 05.70.Ln, 47.70.-n, 87.23. Ge 



I. INTRODUCTION 

Statistical mechanics has revealed many interesting 
properties for nonequilibrium systems. The fluctuation 
theorem (FT) was first discovered for a deterministic 
nonequilibrium system driven by a shearing force [l|-|Jl. 
The FT was also shown to govern a wide class of stochas- 
tic nonequilibrium systems [5|-[y]. The violation of the 
fluctuation and dissipation relation is also an important 
characteristics for nonequilibrium systems J10l4l4| . The 
existence of nonzero current is a main consequence ob- 
served in the nonequilibrium steady state (NESS) closely 
related with the violation of the detailed balance (DB) 
which have been studied up to date |15l - l22| . 

There are several sources to drive a system into 
nonequilibrium such as an external driving 23|, 24 ja con- 
tact with multiple heat or particle reservoirs fl?!, IS, i2(J, 
l2ll . [25|, [2g, a non-conservative force ^ Jj), 22, 27, 2&L a 
time-dependent perturbation for external parameter |29l - 
|35| . etc. We recently studied the case in which nonequi- 
librium is driven by the combination of two sources, the 
nonconscrvativc force and multiple noises that are not 
identically and independently distributed [l^ [22] . There 
may be a stochastic process where thermodynamic quan- 
tities such as energy heat, and work are not defined, most 
likely in biology [3g|- One of us (P. A.) suggested a way 
to transform a stochastic differential equation for posi- 
tion like variables to that for position-momentum pairs 
in zero mass limit 37[. In one dimension it is equivalent 



to the usual overdamped limit while it is not the case in 
high dimensions, giving rise to a new stochastic process 
not derived by the Ito or Stratonovich type calculus [M 
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It is yet an open problem how to construct the prob- 
ability density function (PDF) for the NESS which is 
the counter part of the Boltzmann factor for the equilib- 
rium case. There have been only a few examples where 
the PDF is rigorously found [13, E [13, [M |4l| . The so 
called potential landscape function is defined as minus 
the logarithm of the PDF, which corresponds to inverse 
temperature times the potential or energy for equilib- 
rium systems. Much efforts were made in earlier works 
by Graham et al. [i^j to develop a general formalism to 
construct the potential landscape function, which is not 
complete. The potential landscape with probability flux 
drawn gives a useful and intuitive picture particularly in 
biology, chemistry, and evolutionary science |43| . In this 
paper we obtain the potential landscape function and the 
circulating NESS current for the system driven by a non- 
linear drift force in the presence of noises not identically 
and independently distributed. We flnd a very peculiar 
landscape picture that cannot be expected easily from a 
given drift force. It is the generalization of our previous 
study for the linear drift force [1^ . For the nonlinear drift 
force, which cannot be treated rigorously, an approxima- 
tion method is required. We will use the perturbation 
theory based on the exact result for the linear drift force 
as unperturbed basis, which is valid in the low noise limit. 
The recent theory of large deviation, solving the instan- 
ton solution to cxtrcmize the rate functional, can also be 
applied for the low noise cases [4l|, [4J] . We expect our 
approach to be efficient provided that the unperturbed 
result is known. The two approaches will essentially give 
rise to an equivalent result. 

We organize the paper as follows. In Sec.[n|we discuss 
general aspects on the NESS for the stochastic dynamics 
driven by the nonlinear drift force. In Sec. lIIII we develop 
the first order perturbation theory for the nonlinear drift 
force and obtain the correlation functions by using the 
Fourier transformation of the Langevin equation into the 



frequency space. In Sec. IIVI we find the potential land- 
scape function in compliance with the correlation func- 
tions found. In Scc.|V]we give an explicit expression for 
the NESS current. In Sec. IVII we investigate an example 
for the motion in a double well potential in the presence 
of correlated noises. In Sec. lVIII we present the numerical 
calculation for the perturbation theory and the simula- 
tion for the example. In Sec. IVIIII we summarize our 
main results. 



II. GENERAL ASPECTS 



We consider the Langcvin equation: 
X = f (x) + ^ . 



(1) 



where the state x = {xi, . . . ,Xny , the drift force f = 
{.111 ■ ■ ■ ■, fnY ; and the noise ^ = {ii,...,UY are n- 
dimensional vectors. The superscript t denotes the trans- 
pose of vector and matrix. We consider x to have the 
even parity in time reversal like position variables, ^(t) 
is taken to be Gaussian with zero mean and variance 
{i{tW{t')) = 2Di5(t - t'). D is a n X n real symmet- 
ric matrix, called the diffusion matrix, and is positive 
definite. We consider the noise components ^^'s to be 
not identically and independently distributed. Then D is 
in general not proportional to the unit matrix, which is 
a possible realization of the contact with multiple heat 
reservoirs. We consider the localized system in which f 
has stable fixed points and goes to oo as |x| — > oo. 

The PDF p(x, t) associated with the Langcvin equation 
^ satisfies the Fokker-Planck equation 

^p(x,t) = V.(-f + D.V)p(x,t). (2) 

The steady state is reached when dp/dt = 0. The PDF 
for the steady state reads 



p(x) ex exp (— $(x)) 



(3) 



We call $(x) the potential landscape function. It is equal 
to the inverse temperature times the energy for the equi- 
librium case. 

The probability current reads j — (f -I- D • V^)^. Equi- 
librium is defined as steady state with the DB. The DB 
for even parity variables is found to hold if j = [43 , 
i.e. f = — D • V<I>. In this situation, the drift force is 
exactly cancelled by the diffusive force D • V$ given by 
noises. The NESS is characterized by nonzero j which is 
divergencelcss since p = 0. For the localized system, that 
we are considering in this paper, circulation is the only 
way for divergencelcss current. On the other hand, for 
the extended system it is directed through the system 
from one to another boundary. There were interesting 
works on the noise induced directed current in an one- 
dimensional extended system with the periodic bound- 
ary condition, which explains well the transportation of 
drugs or molecules in biological systems [15|, [la] . There 



arc common ingredients for the two cases, circulating and 
directed currents, which we will discuss later. 

Let us define the force matrix F and the potential ma- 
trix U respectively as 



F,,, 



v,/;, t/,, 



V,Vj$ 



(4) 



Then the DB condition gives U = — D ^F, which leads to 
the condition: 



FD-DF^ = 



(5) 



which can be obtained by using U = U^. 

The thermodynamic equilibrium is given by f = 
— 7Vi? for energy E and D = (^kBT)\ for unit matrix I. 
Since F = F'^, the DB holds. In fact the DB always holds 
for one dimensional systems. For higher dimensions there 
arc two factors that make the DB violated. One is that 
F is asymmetric, F 7^ F'^, which is the case in which f 
is nonconservative, i.e., not derivable from a scalar func- 
tion, f 9^ -W. The other is that D 9^ I. Then Eq. ^ 
may not be satisfied. 

With the broken DB, the system will reach the NESS 
after a long time. Let us define the residual force as 



f,es = f + D • V$ , 



(6) 



which vanishes for the equilibrium case. The NESS cur- 
rent is given as j = iresP- The condition V • j = gives 



V$ • L, 







(7) 



We briefiy summarize our previous study on the case 
with linear drift force, f = F -x, which is known as a high 
dimensional Ornstein-Uhlenbeck process ^M,- The PDF 
is Gaussian and the potential landscape function is given 
by $ = x^ • U • x/2. We found 



= -Q- V$ 



(8) 



where Q is an antisymmetric matrix and independent of 
x. It satisfies the steady state condition ([7]). Then the re- 
sultant current circulates on the cquiprobability surface, 
which can be shown by noting that Q • V<f> is perpen- 
dicular to V$. V$ = U • x yields U = -(D + Q)-^?. 
The condition U = U"^ gives the matrix equation for the 
antisymmetric matrix, 



FQ + QF^ = FD-DF^ 



(9) 



The exact solution was found by using the Jordan trans- 
formation for an arbitrary asymmetric matrix F. The an- 
tisymmetric matrix Q is the single measure for the NESS. 
If Q 7^ 0, the DB is violated, as seen in Eqs. ©, (H]). 
Recently Filligcr and Rcimann [43] showed there exists a 
non- vanishing torque perpendicular to a two dimensional 
heat engine, which is nothing but a manifestation of the 
circulating current in our study. 

For the nonlinear drift force, however, the antisymmet- 
ric matrix is not a sufficient measure for the NESS. We 
can write 



Q • V$ + f, 



off 



(10) 



A new term foff gives a current flowing off tlie equiprob- 
ability surface. Tlien the drift force is decomposed into 
three parts: 



f = -D- V$-Q- V$ + f, 



off 



(11) 



A similar idea of decomposition was used in earlier works 
by Graham and Tel [42, where they concentrated only 
on the dissipative part, the first term, in our terminol- 
ogy. A vector potential A was introduced by Qian [43 
to describe a circular flux and the global transport in a 
two dimensional motor protein movement, where V x A 
corresponds to the circulating current in our study. 

One can see that the probability maximum where 
V<f> = shifts from the fixed point where f = if 
fo// -/L Q n jg contrary to a usual observation that the 
fixed point coincides with the probability maximum, as 
in the DB case where f = — D~^ • V<i>. Q may become de- 
pendent on X. Then the resultant current — (Q(x)- V$)p, 
flowing on the equiprobability surface, may not be diver- 
genceless. However, the total NESS current with ioffP 
added goes divergenceless and circulates off the equiprob- 
ability surface. Interestingly there may appear various 
centers of circulation, which will be seen in the following 
sections. 



III. PERTURBATION THEORY 

Our aim is to find the PDF for the NESS. The nonhn- 
ear drift force cannot be treated rigorously for a general 
D. We use the perturbation theory for the nonlinear 
force expanded about a fixed point. Let x = be the 
fixed point. Then one can write 



ft = F>j'xj + GijkXjXk + HijkiXjXkXi 



fl2) 



The Einstein convention is used for repeated indices de- 
noting the summation. The potential landscape function 
can also be written as 



$ = hiX. 



+ ^ (^^ + u^r) -^-^ 



+ 7y^'^tjkXiXjXk + —A.jkiXiXjXkXi H . (13) 

The linear term in the drift force yields the unper- 
turbed basis for the perturbation theory where the un- 
perturbed potential matrix U'*'^ is known for F^''^ from 
our previous study [19[. The superscript (1) denotes the 
perturbation due to the nonlinear drift coefficients. The 
perturbation theory is expected to be the same as the low 
noise limit, i.e., the small D limit, since the PDF is dom- 
inant near fixed point in this limit. We then proceed the 
perturbation theory up to the first order in Gijk, Hijki, 
both treated as of the same order. It is the same sense as 
in the perturbation theory for the anharmonic effect of 
the lattice vibration where the cubic and quartic terms 
in the potential are treated as of the same order in order 
to give a bound potential. 



Note that the linear term hiXi in $ is responsible for 
the shift of the PDF maximum from x = 0. The moments 
of the PDF comprise a set of correlation functions written 
in terms of the coefficients of $. On the other hand we 
will find the moments directly from the Langevin equa- 
tion in the Fourier space. Then we can determine the 
coefficients of $ in a selfconsistcnt manner by comparing 
the moments found from the two means. 

Let Cij... be the connected correlation function 
{xiXj ••■)£, known as a cumulant. (• • • ) denotes the av- 
erage evaluated by the PDF associated with $. For ex- 
ample, Cij = (xiXj) — {xi){xj). Up to the first order we 
find 



C, 

a 



^il^l (^ i Imn -ti-il -Limn : 

i(i)[ 



(14) 



V — ^ij ^ [f^U R]ij — -AklmnRikRjlRmn , (15) 



a 



ijk 



^ lriin-^il-^J7n^kn 



Cijki — —Aiij>k'i'Rii'Rjj'Rkk'R> 



where we use 



R=(U^"0 



(0)n-1 



(16) 
(17) 



(18) 



One can use the standard Feynman diagram technic for 
vertex tensors f, A to obtain the above results. 

Using the Fourier transformation one can see that the 
Langevin equation, the stochastic differential equation in 
time, turns into the algebraic one for x(cij) in frequency 
space. Up to the first order in Gijk, Hijki we find 

Xi{uj) ^a^j{uj)\^^j{uj) 

/I 
— =La;/'(a;i)a„iTO'(w - uji)^i, {iOi)^rn' {^ ~ ^i) 
V Stt 



27r 



X 6'('^l)Cm'('^2)^n'(w - Wi - W2)| , 

where we introduce a matrix 

a{uj) = \-iuj\ - F(°)] 



(19) 



(20) 



From this we can find the correlation functions 
{xi{ijj)xj{u!') ■ ■ ■) in frequency space, using 



{Uuo)ij{Lo'))=2D,j5{u: + u:') 



(21) 



In this calculation (• • • ) denotes the average over noises. 
The connected correlation functions at equal time can 
then be found by 



C,,...{t) 



diuduj' 
27r 



{x4Lu)x,{uj')---),e-''^'^+^'+-^' . 

(22) 



The connected correlation functions in frequency space 
can be found to have a delta function factor, S{uj + uj' + 
•••). Therefore the connected correlation functions at 



equal time arc independent of time and coincide with the 
corresponding correlation functions in the steady state, 
given in Eqs. (fT4 | -(fT7 l) . 

After a careful calculation, we find the correlation func- 
tions: 



where 



Ci = ~2F- ■ GjlmPlrni^) , 



(23) 



^{a,fe(-aj)/3j,(a.) + (i ^ j)} , (24) 



Ciiu=^Gimn I — — j|Am(w)/3j„(a;')afe/(-w- w') 



dujduj 
72^ 

+ (jofc) + (iofe)} , 

(koduj'duj" r 
7r)3 I 



(25) 



J (27r)3 
(iia{uj)(ijb{i^')l3kc{i^")aw{-uj -J - J') 

+ (fco/) + (jo/) + (^oO} • (26) 

where (i O j) denotes the term obtained by exchanging 
indices i and j. We introduce another matrix 



/3(lj) = a(w)Da^(-cj) . 



(27) 



We can determine the coefficients in $ by comparing 
Eqs. (I23l)-(l26l) with Eqs. (fT4 l) -([T7 | . 



IV. PROBABILITY DENSITY FUNCTION 

Comparing Eqs. (fT5|) and (|24|) . we can find the poten- 
tial matrix in zeroth order 

R^[U("']-i=|^/?(c.). (28) 

This gives the antisymmetric matrix in zeroth order, 

Q("' = -F(") / '^(i^\ + F("))-ii?(zc^/ - f(o)^)-i - D , 

(29) 
where the relation U^") = -(D -H Q(o))-if(°) is used. It 
is very plausible to get the integral representation for 
the antisymmetric matrix which was found as a compli- 
cated series form in our previous work for the linear drift 
case [l9l |. 

For simplicity of notation, let us define 

d(w) = U("'a(cj) , (30) 

/3(c.) = U(°)/3(c.) . (31) 

Comparing Eqs. ([TB| and (P5|) . we find the coefficients 
for the cubic term in $: 



-L ijk — ^ijk '-'abc , 



(32) 



'^abc 
^ijk 



diodiO 



8 / -— T^|/3i6(a;)/3jc(a;')afca(-w-a;') 



+ {j f> fc) + (i o k)\ . 



(33) 



Comparing Eqs. p7|) and (pS)) . we find the coefficients 
for the quartic term in $: 



Ayfc; = *.yfei i?a&cd , 



(34) 



where 



^ ijkl 



dudu'duj" 



{ 



/3ib(a;)4c(w')/3/cd(w")a;a(-i^ - w' - w") 

-h (fc ^ + (j ^ + (* ^ } ■ (35) 

Comparing Eqs. ([15]) and ([M)) in the first order, we 
find the first order correction for the potential matrix: 



r(i) 



1 



where 



^ ij r) \^ijmn i ^ ij J^klvin J ^nin : 



duj 



(36) 



^ = 12/ ^ {hMa,u{-u:) + (* o j)} . (37) 

Finally, comparing Eq. (|14l) and (P5|) . we find the co- 
efficients for the linear term in $: 

K = - ([{D + Q^"))-i],/G,„„ + ir,™,") i?,„„ . (38) 

Then we can estimate the shift 'X-max of the probability 
maximum from x = 0, which is given from V$ = 0, 

\^^nax\i — ^ij j 

-[(F(°))-].G.„„ + ii?.r.™„)i?™„(39) 

From Eqs. ((T4| and (p9)) . we can find (x) as 

(x,)=-[(F("))-i],,G,™„i?,„„. (40) 

It is rather surprising that the shift of the probability 
maximum might not be in the direction of (x) . 



V. NONEQUILIBRIUM STEADY STATE 
CURRENT 



We divide the residual force in Eq. (fTO|) into the zeroth 
and first order as ires = ires + ires- Each term is given 



fW = -Q(") • va> , 

4^1 = C' + fS = -Q''^ • v*'"^ + fo// • (41) 



fori is the term with the first order antisymmetric matrix 
Q'-^-'. toff only appears in the first order, so equals to 



f(i) 



Similarly the NESS current is decomposed as 



;(n) 



i(i) 



J =J' ' +J' ' =J 

= f(0) ^f(l) 
'-resr ' on r 



(0) , !(1) 
' Jon 



Jo// 



(42) 



j*^*^) flows on the equiprobability surface. It is divergence- 
less because Q*^^-' is constant in x. j|,„ also flows on the 
equiprobability surface, while it may not be divergence- 
less since Q'^-* generally depends on x. y^L flows off the 

equiprobability surface. The total current j*^^^ in the first 
order should be divergenceless. The existence of nonzero 
current is the signal for the NESS. The existence of the 
current flowing off the equiprobability surface is resultant 
from the combination of the nonlinearity of the drift force 
and the nonuniformity of noises, which is not seen in the 
linear case. 

The velocity field irel for the total current j^^^ in the 
first order can be found from 

V$(i) = -(D + QW)-i(f(i)-fW) , (43) 

where $^^-' is the first order correction to $ and f'^' is 
the nonlinear part of f . We find 

[f^l)]. = [D + Q(«)],,/^, 
+ [(D + Q(o))U«],fcXfc 



G^kl + -[^ + Q^%jTjkl \xkXl 



Hik 



3! 



[D + Q(°)],yAjfeim 1 XkXlXra ■ (44) 



The resultant current j^^^^ is expected to be divergence- 
less. 

Note that fres in the above equation is independent 
of Q*-^'. In fact there is no unique way to determine 
Q*-^-* and f^j^^ simultaneously. There is a gauge invariance 
under the change: Q(i) ^ Q^^' + SQ^^\ i^^fj -^ iffj: - 
Q(i) . V^^^-*. One way is to extract the antisymmetric 
matrix maximally from ires- We expand the perturbed 
antisymmetric matrix Q^^' as 

Qij + qijkXk + qijkiXkXi -\ , 






(45) 



where gy... :== -qj, 
in Eq. ^t^ reads 



The part to be extracted from f, 



(1) 



(0). 



'qijkiU],^XkXix. 



(46) 



We find 



Qijk = -^ [GikiRi, - {i^ j)) 



Qijkl 



-^ [HikhnRni] - {i ^ j)) 



-— ([D + Q^%,,AnklmR,nj - (« ^ j) 



(48) 



(49) 



Here we use the property that (A — A'^)/2 is the anti- 
symmetric part of a matrix A. In general the resultant 
current yon is not divergenceless. 

The current yL Sowing off the equiprobability surface 

is given from the remaining part in f^es after the maximal 
extraction of Q^^^ We find the velocity field of it as 

[{offh = [D + Q^'%h, 



{(D + QW)UWR + RU(i)(D - Q(o))} U(°) 



Xk 



7:Giki 



rakl 



-i[D + Q^%,nTmkl'Rv^Uf ) XkXl 

T:Hikim + -^H jkim' Rm' iU jj-^ + t;^[D + Q 



(0)1 



\in'—^nkLm 



-^[D + QW],„A„,,„,i?, 



,U. 



(0) 



■in'i^ jjYi ] •^k-^l^m • 



(50) 



•M 



The resultant current yJs may not be divergenceless, 
while the total current of the first order, j^^^ = ji,/ 4- 
JqA, is expected to be divergenceless. As observed in 
an example in the next section, j*-^' shows interesting 
behaviors. In particular the current yields circulation at 
multiple centers which arc not coincident with either the 
probability maximum or the fixed point. 



VI. MOTION IN A TWO-DIMENSIONAL 
POTENTIAL WELL 

Let us consider a motion in a two-dimensional double 
well potential V{x, y) for xi — x and X2 ~ y. The double 
well potential is given by 



y{x,y) = ^^^i^ 



^2 2 



(51) 



There are three fixed points where 'VV ~ 0: (0, 0) (sta- 
ble), (a/2, 0) (saddle), (a, 0) (stable). The diffusion ma- 
trix is taken to have nonzero off-diagonal element D12 . 

We carry out the perturbation theory around (0,0). 
Note that the force matrix, given by Fij = —d'^V/dxidxj, 
is symmetric. However, noise correlation leads the sys- 
tem to nonequilibrium. In zeroth order we get 



9.. = -2 



(D + Q(°))U(i)r1 -(ioj) 



(47) 



F(o) = _ 



ki 

k2 



(52) 



The coefficients of nonlinear terms arc given by 

Htjk = fsSiiSjiSki ,Gijki = /4<^ii (5^1 41(5/1 , (53) 

where 

3fci _ 2fci 

73 — , /4 — 5" ■ 

a a^ 



and 



The full force matrix is given by 



-fci + 2/3X1 +3/4.T?/2 

-fc2 



(54) 



Then the condition for DB, FD — DF'^ = 0, gives 

(-fci + fc2 + 2/3,Ti + SfixDDu = . (55) 

The DB is violated for D12 ^ 0. 

In proceeding the perturbation theory, we need to 
carry out multiple integrals for functions given in terms 
of matrices a(a;), /3(a;) defined in Eq. ([?7)) . We find 







a{io) = ( --+^-- _^ 



— iuj-\-k2 



(56) 



D12 



I {-iu:+ki){iu:+ki) ( — icj+fci )(?:cj+fc2) 

^12 D22 

{ — iuj-\-k2){iuj-\-ki) { — iu-]-k2){iu-\-k2) 



B(uj) = {-■i'^+K:i)(iuj+Ki) (~-iuj+Ki)(t.uj+K2J /-gyv 



One can find U'"' and Q'"-' from the integral representa- 
tions in Eqs. (1^51) and 



(58) 



Dii 


2-D12 


fci 
2D^2 


kl+k2 

D22 


fei + fca 


fe2 



and 



^(0) _ fci - fc2 

^12 - T n— ^12 

Kl + fc2 



(59) 



For simplicity we consider the case where ki = k2 ^ k, 
Dii = D22 = d, D12 ~ ed. Then we find 



Q(0)=0, UW = 



k 



1 -e 



(60) 



d(l-e2) y-e 1 

After carrying out multiple residue integrals, we find 
6 






d(l-e2) ' 



^-"21 — *-'l2 



3e 



2d(l-e2) 



and 



Ti 



11 



2/3 



rf(l-e2) 



r2ii — ri2i — Til 



2e/3 



3d(l - £2) 



(61) 
(62) 

(63) 
(64) 



6/4 



^""- rf(l-62) ' 

A2111 = A1211 = A1121 = A1112 

3e/4 
2d(l-e2) • 

Now U^-'-'' can be found as 

3eV4 






rjil) 



2fc(l-e2) ' 
3e/4 



4fc(l-e2) 



C/2« = 



(65) 
(66) 

(67) 

(68) 
(69) 



The coefficients hi for the linear term in $ can also be 
found as 



hi 
h2 



2^2/3 
3(1-62) 

2e/3 
3fc(l-e2) 



(70) 
(71) 



The shift of probability maximum from the fixed point 
is found as 



-Riihi = , 



Vv 



R2ihi — — 



2edh 
3fc2 



(72) 
(73) 



The shift is made perpendicular to x direction where the 
nonlinear force is applied. It is somewhat contrary to our 
expectation. There may be an effective transverse force 
due to noise correlation. Note that the shift depends on 
D12 and /3. It is a novel phenomenon due the combi- 
nation of noise correlation and nonlinearity in the drift 
force. 

The average position (x) can be found from Eq. (|40)) . 



Note that it is in x direction, perpendicular to yimax- 
The potential landscape function is then given as 

*(^'' y) "= ^^(^' + y' " 2ea;y) + ^{-ex + y) 

ZdK 6K 

-F— — (-ex + xy) + irry-^ + ^^ v) 
4:kK 6aK 



(74) 



(75) 



where k = 1 — e2 Note that the perturbation theory 
is valid for small d where we can treat x to be 0{Vd). 
Therefore the first order perturbation theory with respect 
to nonlinear coefficient f^, /4 is equivalent to the expan- 
sion of $ up to the first order in d. This implies that f^, 
fi are not necessarily small if d is taken to be small. 



The antisymmetric matrix in the first order is found 



as 



A\' 



-|(/3- + /4-^) 



(76) 



As noticed it depends on x. The resuhant current jon = 
— (Q(i) • V<I>*^°^)p flows on the equiprobability surface, but 
is not divergenceless. 

The velocity field for the current j^A flowing off the 
equiprobability surface is found as 



C/i 



rf (1) 1 



-g(.x^-x^y), (77) 

2^^^^+^((l-2e> + .,) 



3fc 
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The velocity field for the total current j^^^ in the first 
order can be found as 



rf(i)i 

I resi 



3edl 



AkK 
3e/4 



-{-ex + yj) 



2e/3. 2 
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(79) 



rf(i 



)1 
es\y 



^e^fs , 3ed/4 /,-, „ 2n , \ 



2e/3 . 2 

(x 



3e/4 



3^ ,^ ex,) - :^ (x^ - .x^y) .(80) 

One can show j^^^ is divergenceless, i.e., V • j'^' = 0, 
as expected from the steady state condition. Therefore 
we expect that it yields circulation off the equiprobability 
surface. The current vanishes at the center of circulation. 
Surprisingly there are two centers at (xc,yc)' 



Vc 



(81) 



where we assume d is small, which is the criterion for the 
perturbation theory. 

VII. NUMERICAL STUDIES 

In order to confirm the results obtained from the per- 
turbation theory we solve the stochastic differential equa- 
tion ^ numerically. We consider the example in the last 
section. Let {t„;n — 0,...,N} be discrete time steps 
with interval At. We write x(i„) = x„, 2/(t„) ~ yn. 
Then x„ and j/„ are updated as 



Xn = x„_i + /^,„-iAt + (^Vl - £^6 + e6) V2dAt , 

yn=yn-l+fy,n-lAt + ^2V2dAt, (82) 



where f^^n-i = -kxn^i + f3x'i_i + /4a;^_i, fv,n-i 
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FIG. 1. The contour plot of the PDF from the numerical 
simulation for d = 1.0, e = —0.3, a = 2.0, k = 1.0. 



from {—1, 1} at each time step. We execute the simula- 
tion up to A^ = 10^ steps. 

The shift of the PDF maximum is proportional to d, 
as shown in Eq. ([73|). In order to show the shift clearly 
we take d = 1.0, that is beyond the perturbation theory. 
Fig.[T]shows the contour plot for the PDF obtained from 
the simulation. We use a negative value of D12, so the 
shift from (0, 0) is in +y direction, agreeing with Eq. ([75]) . 
Around the saddle point (1,0), where /m = 0, there is 
no shift, which can be seen from Eq. (pQ]). The shift 
from the other fixed point (2, 0) is in opposite direction, 
as can be expected from the symmetry of the potential 
well. Fig. [5] shows the two contour plots of the PDF ob- 
tained from the perturbation theory and the numerical 
simulation. We take d = 0.1 for which the perturbation 
theory is relevant. The two plots seem to be well coinci- 
dent to each other near (0,0), maybe within the window 
of side ^ Vd, that is about 0.3. In Fig. |3] we present 
the PDF along the two lines, a; = and y = 0, through 
the origin. The maximum of P(0, y) is the shifted value 
of the PDF maximum. It is estimated from the pertur- 
bation theory as 0.05. The solid line is the plot from the 
theory and the scattered circles from the simulation. The 
figure shows a good agreement between the two methods. 

Fig. |4] shows the NESS current drawn on the con- 
tour lines. For this example the current in the zeroth 
order vanishes. We take a larger value for a. The 
larger a, the smaller /a and f^. Then the pertur- 
bation theory will give a more accurate estimation in 
the region far from the origin. The figure on the top 
left shows jon = — (Q^^-* • V<I''"-')p flowing on the con- 
tour line (equiprobability surface), the figure on the top 
(1) 



^1^2 are random numbers chosen independently right shows j^ A = ioffP flowing off the contour line. 




Current on the contour 



Current off the contour 



FIG. 2. Contour plots for PDF from theory and simulation 
for d = 1.0, e = -0.5, a = 2.0, k = 1.0. 
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FIG. 3. PDF along a- = and y — lines through the origin 
for the same values of parameters as in Fig. [2] The maximum 
of P(0, y) is estimated to be 0.05, which is the value of the 
shift of the PDF maximum from the origin. The solid line is 
for the theory and the red circles for the simulation 



and the figure at the bottom shovifs the total current 
j'-'^' — 2o7i +Jo/f The current on the contour line shows 
a cut (discontinuity line), at cc = 0; the current flows into 
the cut for y > and flows out from the cut for y < 0. 
It manifests that the current is not divergenceless, as ex- 
pected from the dependence of Q^^^ on x. The current 
jLf is not divergenceless, not clearly seen in the figure. 
It can be observed that there is a local flux flowing down 

through the origin, filling up the discontinuity of jon • As 
a result, the total current becomes divergenceless while 
circulation splits into two parts. It is not an artifact of 
the perturbation theory, for the PDF is accurate in the 
region near the centers of circulation. For d = 0.1 and 
k = 1.0, they are located at (±0.316, =fO. 158), estimated 
from Eq. (|81l) . The circulation takes place across the 
contour line, which is an important characteristics for 
the NESS together with the shift of the PDF maximum 




FIG. 4. NESS Current. The current flowing on the contour 
(top left), the current off the contour (top right), the total 
current (bottom) . We take the same values of parameters in 
the last figure except for a = 10.0. There appear two centers 
of circulation for the total current at (±0.316, =f0.158) 



from the fixed point. However, it is quite surprising that 
there appear several centers of circulation. 



VIII. DISCUSSION 

We have investigated the potential landscape and the 
NESS current for the localized stochastic system driven 
by the nonlinear drift force and in contact with nonuni- 
form noises. Unlike the extended system the induced 
NESS current is circulating. For the linear drift force 
it was found to circulate around the equiprobability sur- 
face [1^ with the probability maximum at the fixed point 
of the drift force. In this study we have found that the 
NESS current fiows off the equiprobability surface and 
circulates around multiple centers. The PDF maximum 
is found to shift from the fixed point of the drift force. 
For an example we take a familiar potential well yielding 
a conservative force. Difference is made only by noise cor- 
relation. However, it turns out that the resultant poten- 
tial landscape, having multiple centers of current circula- 
tion, is rather drastically deformed from the equilibrium 
shape. Circulation is found to split into two parts and the 
boundary current fills up the discontinuity of the current 
flowing on the equiprobability surface. However, further 
investigation is needed for a more plausible explanation 
based on the physical origin. Noise induced transverse 
force might be a possible origin, which is not confirmed 
yet. They are novel phenomena caused by the combina- 
tion of two causes, the nonlinearity of the drift force and 
the nonuniformity of high dimensional noises. Interest- 
ingly they arc basically equivalent to the two ingredients 
for the noise induced transportation current found for 
the extended periodic system: the ratchet type potential 



(nonlinear force) and the additional noise (high dimen- 
sional nonuniform noises) jlSl |l6[ . Our finding is based 
on the first order perturbation theory that is valid in the 
low noise limit. It is confirmed by numerical simulations. 

It is interesting to realize the NESS with a circulat- 
ing current experimentally, probably in an optical trap 
experiment in two dimensions. It is a challenging task 
to produce noises experimentally which are not identi- 
cally and independently distributed. Force in experi- 
ments may be conservative and then the force matrix 
F is symmetric. In this case one can still produce a 
nonequilibrium situation if the principal axes of the dif- 
fusion matrix are chosen to be different from those of F, 
as suggested by Filliger ct al. [43|. 

The fluctuation of nonequilibrium work production has 
recently been studied for the linear case [23, Um ■ There 
are interesting properties found such as the exponential 
tail with a power law prefactor of the work distribution 
function and the dynamic phase transition in the expo- 



nent of the prefactor. The properties of the NESS for the 
nonlinear case are found to be quite different from that 
for the linear case. In this sense it will be interesting to 
study the nonlinear case from the point of view of the 
fluctuation theorem. 
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